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The expansion of the universe causes spacetime curvature, distinguishing between distances mea- 
sured along and transverse to the line of sight. The ratio of these distances, e.g. the cosmic shear 
distortion of a sphere defined by observations of large scale structure as suggested by Alcock & 
Paczyiiski, provides a method for exploring the expansion as a function of redshift. The theoreti- 
cal sensitivity to cosmological parameters, including the dark energy equation of state, is presented. 
Remarkably, sensitivity to the time variation of the dark energy equation of state is best achieved by 
observations at redshifts z < 1. While systematic errors greatly degrade the theoretical sensitivity, 
this probe may still offer useful parameter estimation, especially in complementarity with a distance 
measure like the Type la supernova method implemented by SNAP. Possible future observations 
of the Alcock-Paczyriski distortion by the KAOS project on a 8 meter ground based telescope are 
considered. 



I. INTRODUCTION 

We now have strong evidence that the expansion of the 
universe is accelerating, from the original method of Type 
la supernova distance-redshift measurements |^ and 
concordant observations of the cosmic microwave back- 
ground (CMB) power spectrum and of large scale struc- 
ture H . Understanding the nature of the dark energy 
responsible for the acceleration will have profound im- 
plications for cosmology, high energy physics, and fun- 
damental physics. Mapping the expansion history of the 
universe offers a way to gain insights into the dark en- 
ergy and the fate of the universe, for example by charac- 
terizing the equation of state behavior which is directly 
related to properties of the scalar field potential. 

Distance measures, notably the supernova method, 
have proved useful at constraining the energy density and 
equation of state of the dark energy, with great improve- 
ments expected in the next decade. But these involve 
an integration over the expansion rate behavior H{z), 
which itself involves a redshift integral over the equation 
of state w{z). We can ask whether we can devise a more 
direct probe of the acceleration. In fact one such, the 
redshift drift test, was proposed by Sandage [D in 1961 
and developed further by Linder |6). 

The redshift of a source is a central astrophysical ob- 
servable. It is directly related to the change in time in- 
tervals due to the cosmic expansion between a photon's 
emission and observation, z = dto/dtf, — 1 = ao/a-e — 1, 
where a{t) is the scale factor of the universe. But obvi- 
ously one could consider a second derivative term, a time 
dependence of the redshift itself as the universe ages: 
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This provides a direct measure of acceleration, being 
effectively a second time derivative, as can be seen from 
the low redshift limit: Az/z « —qoHoAt, where qo is the 
present deceleration parameter. However since the astro- 
nomical observing time is much smaller than the Hubble 



time, At <C H^^, this is not a practical probe. For ex- 
ample in the most optimistic case of observing over a 
period of 10 years a hypothetical spectral line emitted at 
the CMB last scattering surface at z = 10^, one requires 
a redshift measurement of a part in 10^ to distinguish 
cosmological models. For emission line objects at z = 5, 
this becomes a part in 10^. 

Moreover, just as peculiar velocities affect redshift 
measurements at the level of v/c, so do peculiar accel- 
erations, i.e. local gravitational potentials $ from inho- 
mogeneously distributed matter, affect the redshift drift 
measurements at a level ^/c^ ^ 10^^. This latter effect 
even ruins the generalization of the redshift drift called 
the cosmic pulsar test, where timing is improved by mea- 
suring a large number A'^ of wavelengths or pulses . 

But if we are stymied in measuring the acceleration 
directly, at least we can hope to measure the first deriva- 
tive of the expansion, H{z). One of the ways to do this 
is the cosmic shear, or Alcock-Paczyhski test. Section 
H sets up the formalism while Section [II applies Fisher 
matrix analysis to investigate the theoretical sensitivity 
of this method for estimating the cosmological parame- 
ters. In Section IV we introduce observational sanity in 
the form of systematic uncertainties and discuss the pro- 
posed KAOS project as a means of carrying out this test. 
We summarize our conclusions and plans for future work 
m Section 0. 



II. COSMIC SHEAR TEST 

Proper distances measured along the line of sight carry 
information through light emitted at different times in 
the expansion history of the universe. Therefore in a 
sphere of comoving points, differences in the emission 
times lead to probing the geometry at different expan- 
sion rates. So observationally a sphere will appear to 
be distorted, or sheared, with the magnitude of the ef- 
fect sensitive to the expansion rate. This is the cosmic 
shear effect discussed by Alcock & Paczyhski ( ; do not 
confuse this with shear from weak gravitational lensing - 
due to inhomogeneities rather than the global structure 
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of spacetime) . 

In more detail, if we consider the small difference in 
radial distance between nearby emitters, then we local- 
ize the behavior in redshift and essentially measure the 
expansion rate at that time. However transverse to the 
radial direction the angular separation between comoving 
points is measured at a constant value of the scale fac- 
tor a(te), and then the light from each source propagates 
over the same intervening distance to the observer at to- 
Thus in the radial case the data gives a snapshot of the 
expansion rate while the transverse distance contains an 
integration over the expansion history from emission to 
observation. 

Consider a sphere of comoving points. The distance 
through the center of the sphere along the line of sight is 
the proper distance. 



(ir|| ^ dt = a(t) dvc 



The transverse distance is 



{1 + z)-^H-\z)Az. 



(3) 



dr^ =r„(z) A0= (l + z)-V,(z) A0, (4) 

where z = a^^ — lis the redshift, H — d/a is the Hubble 
parameter, Tc the comoving distance, and the angular 
distance (see, for example, |^). 

From the observables of the angular scale of such 
comoving sources, their central redshift z, and their red- 
shift extent Az, one can form a quantity 

D{z) = Az/Ae = H{z) ( dz'/H{z'). (5) 
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The cosmic shear is then 
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This has some excellent properties for a cosniological 
probe. In particular, it has dependence on H{z) directly, 
not just through an integral. Since H is related to the 
energy density of the universe then one can try to map 
the density history and equation of state. Indeed for 
time varying equation of state of the dark energy w{z), 
H depends on an integral of w{z), so distance measures 
involve w{z) as a double integral. Therefore one might 
hope that the Alcock-Paczyhski differential distance test 
might be more sensitive to reconstructing the dark energy 
equation of state than a standard distance test. 

Another interesting characteristic of is that it does 
not depend on any absolute scale, since H appears in 
both numerator and denominator. So there is no absolute 
measurement to marginalize over. 

The physics of the test seems clean, using pure ge- 
ometry of the background spacetime, so long as we can 
find sources that are comoving and defining a known local 
spatial geometry. Conventionally the local source geome- 
try is taken to be spherical, as it would be for an isotropic 



arrangement. Possibilities for defining isotropically ar- 
ranged sources include large coherent objects such as su- 
perclusters in the linear density regime and so following 
the isotropic expansion of the universe, or lengths defined 
through correlation functions between individual objects, 
such as for galaxies or Lyman alpha forest absorbers. 
The extent to which these assumptions break down or 
cannot be corrected yield systematic errors to plague the 
method. As for many other cosmological probes, the sys- 
tematic errors turn out to be more severe than statistical 
errors from insufficient data. Section IV discusses this 
further. 



III. SENSITIVITY TO COSMOLOGICAL 
PARAMETERS 



We begin by considering a purely theoretical analysis 
of the capabilities of the test, leaving observational real- 
ities for § tV , A good way to understand the sensitivities 



and degeneracies of cosmological probes is through Fisher 
matrix analysis of the dependence of the observable on 
the parameters. We take a flat universe defined by three 
parameters: the dimensionless matter density r2„j (so the 
dark energy density is 1 — i^m), the value of the dark en- 
ergy equation of state today wq, and a measure of its time 
variation w' . Since we want to consider observations ex- 
tending to redshifts greater than unity, e.g. Lyman alpha 
observations are most plentiful with z « 3, we adopt the 
equation of state parametrization w(z) — wa+Waz/{l+z) 
with the definition w' = dw/dhi{l + z)\z=i ~ Wa/2 
This approximates well the behavior of several classes 
of dark energy models, especially those with a slow roll 
phase, is well behaved even for z > 1, and allows insight 
into the effects of the physically expected time variation 
in the equation of state. 

The distance distortion D{z) = H{z) J dz/H differs in 
its behavior in an interesting way from distance or vol- 
ume measures: the two factors actually depend on the 
cosmological parameters in opposite ways because they 
have reciprocal dependence on the expansion rate. At 
low redshift the direct H{z) dependence dominates since 
all distances must be similar. But at high redshifts the 
universe was matter dominated and so H{z ^ 1) is in- 
sensitive to the equation of state parameters and the dis- 
tance factor takes over, since it retains memory of those 
parameters due to its integral nature. 

Figures and |2| show the dependence of the Hubble pa- 
rameter and the distances Tc and r^, respectively, on the 
cosmological parameters fim and wq {wa behaves sim- 
ilarly). One clearly sees that the dependencies are in 
inverse relation. This implies that there can be crossover 
ranges in redshift where the distance distortion D is es- 
sentially independent of one of the parameters. While 
this makes it impossible to determine that aspect of cos- 
mology with observations in that redshift range, it has 
a benefit as well. From the figures one can deduce that 
degeneracies exist where one parameter can be adjusted 
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to counteract the effect of another. But at the crossover 
points one parameter wiU not affect the distortion and 
so the degeneracy can be broken. Essentially, observa- 
tions near a crossover apply to a reduced phase space, 
and hence the other parameter estimates will be sharper. 
This is clearly seen in later figures. 




z 



FIG. 1: The Hubble parameter H[z) as a function of redshift, 
for the fiducial flat model Q.m = 0.3, w = —1 (solid curve) 
and variants. Upper (lower) dotted curves have Q.m ~ 0.35 
(0.25); upper (lower) dashed curves have w — —0.8 ( — f .2). 



An excellent way to explore the sensitivity of an obser- 
vational test to the cosmology is through Fisher matrix 
analysis This methodology approximates the likeli- 
hood surface of the parameter fit to observations with a 
gaussian probability near the best fit, the fiducial model. 
The sensitivity of its estimation of parameter values de- 
pends on the derivatives of the observational quantity 
with respect to the parameter, dD/dx, and the precision 
with which the observations can be made. Writing the 
errors as SD = {SD/D)D, we see that for a given frac- 
tional measurement precision ao = 5D/ D we can obtain 
a parameter error estimate 5x by investigating din D/dx. 
Figure Isl shows this central quantity of Fisher analysis. 
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FIG. 2: The comoving distance Tc. and angular distance Ta as 
a function of redshift, for the same models as Fig. |^. However 
now upper and lower models are interchanged, e.g. the upper 
dotted curve has Q.m ~ 0.25 etc. 



FIG. 3: The sensitivity of the observable, here D = 
H{z) rc{z), to the cosmological parameters is encoded in the 
derivatives plotted here. The larger the absolute magnitude of 
the derivative at a particular redshift, the more constraining 
the observations there. 



Indeed we can read off almost all of our results 
from this sensitivity graph. The previously discussed 
crossovers are clear, with the probe aloof to wq around 
z « 1.3 and to w' (i.e. Wa) around z « 2.3. Sensitivity 
to Qrn grows out to z « 1, after which it levels off; wo 
has a sweet spot at 2; ~ 0.4 and again above z « 2; the 
time variation w' has a broad, though low, impact from 
z Ri 0.6 — 1.1. However at z w 1.3 we expect w' to be 
uncorrelated with wq and so possibly easier to determine. 

Also note that while at low redshift the sensitivities to 
flm and Wq enter with the same sign, this changes after 
the Wo crossover. This means that at low redshift the 
error ellipse in the wq — plane will have the same 
orientation as for the supernova distance case: making 
Wq larger (less negative) can be counteracted by mak- 
ing rim smaller, defining a degeneracy direction in that 
plane. However higher redshift cosmic shear observations 
will have an orthogonal degeneracy direction, holding out 
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the promise of complementarity with supernovae. A sim- 
ilar rotation of error ellipse contours with redshift can be 
predicted between wq and w' . 

Figure ^ for the Fisher sensitivity even allows us to 
calculate lower limits on the precision with which we can 
estimate the parameters: 

Sx> {dlnD/dx)-^(jD- (7) 

This lower limit refers to fixing all parameters but one, 
and so will underestimate the true error due to degenera- 
cies. If we take, say, 1% precision in observations, then 
the lower limits on estimating {ilrn,uio,w'} are 0.015, 
0.1, 0.15 for the peaks of the sensitivity curves. However 
multiple observations can improve on the precision, while 
systematic errors will put a floor on it. To proceed fur- 
ther quantitatively we must input an observational suite 
into the Fisher method. 

To best illustrate the results we begin with a simple set 
of observations. We assume equal numbers n of obser- 
vations in redshift bins of width 0.1, each with the same 
precision ct_d, and vary the redshift range the observa- 
tions cover. Initially we apply only statistical errors and 
so only the combination ao/y/n matters. 

Purely to test our intuition from the Fisher sensitivity 
figure, we take the highly idealized situation of constant 
w{z) (i.e. fix w' = a priori) in Figure ^ and fixed Qm in 
Figure |^, as well as unrealistically good precision. The 
rotation of the contours with redshift is clear. Note that 
for the crossover redshift z w 1.3 the determination of 
wq is uncorrelated with the other parameters, i.e. the 
contours are vertical or horizontal. 



IV. SYSTEMATIC UNCERTAINTIES AND 
OBSERVATIONAL METHODS 

Virtually all proposed observational methods for prob- 
ing cosmology run into limits imposed by systematic un- 
certainties rather than statistical errors. So although the 
cosmic shear test appears rather promising in its sen- 
sitivity to parameters and complementarity with other 
probes, as expected from its inclusion of a bare factor 
H{z) as discussed in the first two sections, it behooves us 
not to make estimates of its power merely by speculating 
on achievable future survey statistics. For more realistic 
assessment of the promise of this method for determi- 
nation of cosmological parameters, we must investigate 
the effect of systematics. Detailed characterization of ir- 
reducible uncertainities requires a comprehensive survey 
design and analysis; instead we present here a simple 
model that should illustrate the main effects and give 
reasonably realistic quantitative results. 

We adopt a precision of 2% in measurement, with 
n = 10 observations per 0.1 redshift bin, but also include 
an irreducible systematics floor of 2% in a bin. This 
makes the actual numbers used for the statistical error 
moot, except when we later consider a systematic that 
declines at low redshifts. Generally we adopt a gaussian 




FIG. 4: Idealized results for estimation of w and Qm, with 
only statistical errors. Note the equation of state is fixed 
to be constant a priori. This is purely illustrative, showing 
the rotation of degeneracy directions and decorrelation at the 
crossover redshift. The sizes of the ellipses are idealized, cor- 
responding to pure statistics for 40 observations of precision 
1% in a 0.1 bin in redshift. 




FIG. 5: Idealized results for estimating w' and wo, with fim 
fixed and only statistical errors. This is purely illustrative, 
showing the rotation of degeneracy directions and decorrela- 
tion at the crossover redshifts. The sizes of the ellipses are 
idealized, corresponding to pure statistics for 40 observations 
of precision 1% in a 0.1 bin in redshift. 



prior on the matter density of 0.03, but also investigated 
0.01. The plots show results for observations over red- 
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shift ranges, e.g. z = 1 — 2. 

Figure ^ shows the w — plane, disaUowing the pos- 
sibihty of any time variation in the equation of state - an 
a priori assumption without justification. Note that sys- 
tematics have a large effect on the Alcock-Paczynski test. 
In particular they wipe out most of the complementar- 
ity at high redshift with the supernova method that was 
given by the w sensitivity crossover and resulting rotation 
of the error ellipse. (Note that the contours for the next 
generation supernova survey, SNAP, include systemat- 
ics.) Some complementarity is retained for observations 
at z > 2, allowing an improvement in determining w and 
^Im by a factor two. Of course modeling the dark en- 
ergy as a constant equation of state becomes even more 
suspect as one increases the range of observation. 
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FIG. 6: Parameter estimation in the w — Qm plane, with 
systematics. Note the equation of state is assumed constant 
a priori. Contours correspond to the 68% confidence level. 

When one removes the constraint that the equation of 
state must be constant a priori, the Alcock-Paczyhski test 
does not offer any improvement to SNAP in estimation 
of the present value of the equation of state wq, as seen in 
Figure ^. The combination of the two experiments does 
help to limit flm, but so does, for example, the weak 
gravitational lensing survey that is an integral part of 
the SNAP mission. 

In the w' — w plane (recall w' = Wa/'2), Figure || illus- 
trates that not fixing has a drastic effect on the error 
contours in comparison to the thin ellipses of Figure |^. 
Uncertainty in Qm broadens the contours in a more or 
less fixed direction in this plane, though, so it affects the 
ellipses for some redshift ranges less than others. For ex- 
ample data around z « 0.7 does not suffer as much loss 
of sensitivity to the equation of state parameters. Indeed 
observations covering the range of z = — 1 at 2% pre- 
cision permit determination of the time variation of the 
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FIG. 7: Parameter estimation in the luo — plane, with 
systematics, marginalizing over w' . Removing the prior that 
w = constant strongly inflates the error contours. 



equation of state with errors about 30% looser than from 
SNAP. Unfortunately the respective error regions mostly 
overlap, with little complementarity. Still, they would 
provide independent checks of the results using very dif- 
ferent methods. 
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FIG. 8: Parameter estimation in the w' — plane, with 
systematics, marginalizing over Q.m- Realistic assessment of 
systematics is key to evaluating the impact of the cosmic shear 
test. 



In a bid to optimize the parameter determination we 
can examine how best to concentrate the observations. 
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and over what redshift range. Given a systematics floor, 
an increase in the numbers of observations within a bin or 
an improvement in the statistical precision accomphshes 
httle. An exception to this would be if external system- 
atics needed to be reduced, for example by dividing the 
data into subsets for a "like to like" comparison. As for 
the redshift range, the results found here indicate that 
relatively little leverage is gained, for the given prior and 
systematics, outside the range z ~ 0.2 — 0.9. 

This is well suited for ground based observations by 
large telescopes within a decade. Such a cosmic shear 
survey could be carried out by the KAOS project: the 
Kilo- Aperture Optical Spectrometer proposed as a front 
end for the 8 meter Gemini South telescope. This would 
have multiplexing capability from some 4000 apertures to 
measure detailed velocity maps of supercluster environ- 
ments. In such a linear overdensity region one might hope 
to apply the Alcock-Paczyhski test without any complica- 
tions of nonlinear gravitational physics or gas dynamics. 
Another possibility is studying the anisotropy of corre- 
lation functions of subclasses of bright galaxies. Both 
approaches can cover the preferred redshift range and 
statistics would not be a problem with a 1.5 square de- 
gree field of view and coverage of some 400 square degrees 
of sky to measure precise redshifts for 10^ galaxies. 

Observations at z > 1 were shown to be fairly insensi- 
tive to the dark energy equation of state for the cosmic 
shear test, and thus can be used robustly to learn about 
astrophysics, perhaps from Lyman alpha forest observa- 
tions. Of course the rich panoply of data from such a 
next generation survey as KAOS can be examined with 
other cosmological probes as well (with similar cautions 
and care for the influence of systematic uncertainties). 

V. CONCLUSION 

The cosmic shear test looks extremely promising the- 
oretically for the determination of cosmological parame- 



ters. This is evident from its tomographic dependence on 
the expansion rate of the universe, shown by the appear- 
ance of the Hubble parameter H(z) by itself. It has fur- 
ther interesting properties in the redshift evolution of its 
parameter degeneracies and complementarity with other 
probes. 

However, systematic uncertainties, most probably in- 
volving peculiar velocities and distortions related to the 
local environment rather than cosmic expansion, put se- 
vere limits on the probe's ability to fulfill its potential. 

Unless systematic uncertainties can be brought under 
the 2% level, the cosmic shear test applied at low or high 
redshift does not appear to offer significantly complemen- 
tary or generally comparable limits to the supernova dis- 
tance method. Considering a linear systematic of 0.02z, 
a matter density prior a{flm) = 0.01, or the addition of 
Planck CMB data does not greatly affect these conclu- 
sions. 

The best hope for applying this method does, however, 
lie at redshifts z < 1, which is observationally feasible. 
Therefore careful study of the systematic uncertainties 
might yield some regime in which the next generation 
of wide, deep redshift surveys can bring this method of 
determining cosmological parameters to fruition. 
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